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I have used a newtonian infinite body problem to model a protoplanitary- 
like hot accretion disk and added terms of laminar and Stokes friction. 
Then I used qualitative methods to show the unstability of local regions 
leading to a lemma of radial migration and local unstability in hot, proto- 
planetary-like accretion disks. Then I considered a general infinite body 
problem with a position dependent, gravitation-like force and a general 
perturbating term. 



Introduction 



In (|Hartungl f2009)) I derived an oo-body model of a protoplanitary-like disk mainly 



considering the barycenter of the disk itself and a spherical local region, the observed 
object is located in, with its radius and barycenter. The Lagrangian of a dense region 
in this model is 

2(T-U) _ . 2 , . 2 , 5i , , ({x~x) 2 + {y-^) 



c = -4: — L = * +y + — + "2 ( ^ — ^ ^ 2 - 3 ) (!) 



m 3 



(x 2 + y 2 ) '■ 



with 7713 being the mass of the object, x and y the coordinates with respect to the 
barycenter of the disk, \ an d ip the coordinates of the local barycenter with respect to 
the barycenter of the disk, R the radius of the local region and a; giving the strength of 
the accreting force. This model h as shown to give an explanation to the observations 
described bv lRoskar et al. I ()2008h . 



Changes to the model 



For now it will be interesting to see how the results will change using friction. 
Therefore I will consider laminar and Stokes friction. As for laminar friction the 
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equation of motion reads 



and for Stokes friction it reads 
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with ki s > 0. Thus the systems are 
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Stability of local trajectories 

To retrieve informations on the stability I will search for Lyapunovfunctions Vi. s of 
both systems. They will have to integrate the pde's 



0>Vi := (VV l \f l )=v x d x V l + 
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(x 2 +y 2 ) 2 J 



+ VydyVl + 
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+ a 2 (y - V>) - KiVy d Vy Vi (6) 



and 
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> V a := (W s |/ S ) = v x d x V s + \ +a 2 (x- X )- n s v 2 x d Vi V s 
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+ a 2 {y-^)-K s vl]d Vy V s (7) 
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For Stokes it's easy to see, that 



V s := { 
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with x — x 



x' 2 - v' 2 + y' 2 + v' 2 ) - 

y-ip 



x' > 0, v' x > 0, y' > 0, i£ > 

> 0, v' x > 0, y' > 0, v' y < 

a:' > 0, < > 0, y' < 0, uj, > 

a;' > 0, v' x > 0, y' < 0, v' y < 

x' >Q, v' x < 0, y' > 0, v' y > 

a;' > 0, v' x < 0, y' > 0, ^ < 

x' >Q, v' x < 0, y' < 0, v' y > 

x' > o, < < o, y' < o, < 

x' < 0, v' x > 0, y' > 0, ^ > 

x' < 0, v' x > 0, y' > 0, ^ < 

x' < 0, v' x > 0, y' < 0, ^ > 

x' < 0, v' x > 0, y' < 0, < 

x' <0, v' x < 0, y' > 0, v'y > 

a;' < 0, v' x < 0, y' > 0, ^ < 

x' <0, v' x < 0, y' < 0, v'y > 

x' <Q, v' x < 0, y' < 0, v'y < 
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: also integrates ([7]) on 
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with T4(X; "xi ^vO = = V s (x,v x ,ip,v^) and coordinates respectively chosen to 
compl y x,y,v x , v„,x, i /> > 0, since k s w 2 > < k s v 2 and thus y s < V < (using 



from (jHartung I 



Reducing 17 s to 



:= <j (x,v x ,y,v y ) T ; \x'\ < 
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Vi = V s still holds ©. Obviously £1/ as well as are depending on the coordinates 
chosen. In case the cho sen system gives fli iS = simply rotate the system. This will 
always be possible (see lHartung I (|2009)). In any given system one will be able to 
estimate the global constrained supremum of R holding 

B{{x,^,v x ,v^,x,^,v x ,v^,);R) n f\ s = B((x,tp,v x ,v il) ,x,tp,v x ,v i ,); R) (11) 

B(a;r) be the open ball with center a and radius r. B({x, ip, v x , ty,, x, ipi v x-> v 4>)'i R) 
then is the largest possible local region to be considered. 
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Summerization 



Recapitulating the calculations above and ( Hartung ( 2009h ) they can be summer- 
ized as a 



Lemma 1 

Let X := (x,y,v x ,Vy,x,^,v x ,v^) G X, with 

X := {X; x,y,v x ,v y ,x,ip > 0, 

arg(a; + iy) - axg(x + iy ) = arg(x + iip) - arg(xo + #o) = *P G R} 

and Xq := (xq, yo, v xo , xoi ipo> v xo > u i/>o) the coordinates in any given system of 
coordinates, be the center-of-mass coordinates of a dense object/region D (x, y, v x ,v y ) 
and its local surrounding region L (x, ifr, v x ,v^) with respect to the barycenter of a hot, 
protoplanetary-like accretion disk A. Let D be in the middle region of A. Furthermore 
let a\ be the effective gravitational constant of the center-of-mass of A, let a<i be the 
effective gravitational constant of the center-of-mass of L and let k be the effective 
dissipative, velocity proportional force parameter, then D will radially migrate through 
A and A is unstable on resolution lengthscales R> r holding 

r := sup{e > 0; 3^ G K : B((L, L); e) n % = B((L, L);e)} 



with 
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Generalization 



For consideration of more general scenarios, I will still use the co-body problem 
with the same definition of phasespace-coordinates. But the system shall be altered 
to three dimensions in space and considering a general gravitation-like force F = 
(F x , F y , F z ) T (x, y, z) and a perturbating term 

G = (G x ,Gy,G z ) T (x,y,z,v x ,v v ,v z ,x,'>P,C-i v xi v i>i v <)- Tnus the system reads 
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A Lyapunovfunction V needs to hold 



v z 

F x + G x 

Fy + Gy 

V F z +G z J 



(12) 



> V := v x d x V + VydyV + v z d z V 

+ (F x + G x )8 Vx V+ (F v + G y )d Vy V+ (F z + G z )d Vz V (13) 
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Defining primed coordinates analogous to a;' := x — x it is possible to define V like 

V := cxp (- sgn(x'v x )x' 2 - sgn(y'v y )y' 2 - sgn(z' 'v z ) z' 2 - sgn(v x {F x + G x ))v' 2 

~ S gn(v' y (F v + G y ))v' 2 - sgn(v' z (F z + G z ))v' 2 ) ~ 1 (14) 

Since coordinate functions can be chosen to just be zero at discrete times and finitely 
often on open timeintervals, they do not have any impact on the stability. Furthermore 
accreting systems need to hold (F + G) < at most times. (F + G) > gives the 
same stability implications as (F + G) < 0. Hence it is only interesting to look at 
(F + G) = 0. Suppose this is true for a discrete time then it doesn't have any impact 
on physics like the coordinate roots. As for timeintervals with (F + G) = obviously 
D behaves like a free particle and thus L will not be stable. Hence V holds (fT3|) as well 
as V(L, L) — V(L, L) — with L :— (x, "0, Cj "x,"^,, ), after continuously completing 

V at non-differentiable points. Thus every region around (L, L) contains points with 

V < making the local region unstable and assuring this region to have a radius 
i?>0. 

Implication 

As for any local region L to be unstable implies that any dense object or region 
D will leave L in a finite time. iV-body systems, with N > 2, have trajectories on 
which energy efficient radial migration is possible (for calculations see iDellnitz et al. 
(2006)). It is well known that L migrates on such trajectories. Hence D also migrates 
on such trajectories, but soon leaving L and entering a different local region, which 
will also travel along such trajectories. Thus D has a high potential of migrating large 
distances radially through the accreting system A. The average distance grows with 
R. This will result in a high mixture of objects within A. 

Theorem 1 

Let D be a dense region or object in an accreting system A with local region L, then L 
is unstable, and D will be migrating radially through A. Phenomenologically objects 
D in A will be close to objects D' born far away in L' ^ L. 
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